Molecular-Brownian Correspondence in the Hard-sphere Dynamic Universahty Class 
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We perform systematic simulation experiments on model systems with soft-sphere repulsive in- 
teractions to test the predicted dynamic equivalence between soft-sphere liquids with similar static 
structure. For this we compare the simulated dynamics (mean squared displacement, intermediate 
scattering function, a-relaxation time, etc.) of different soft-sphere systems, between them and with 
the hard-sphere liquid. We then show that the referred dynamic equivalence does not depend on 
the (Newtonian or Brownian) nature of the microscopic laws of motion of the constituent particles, 
and hence, applies independently to colloidal and to atomic simple liquids. In addition, we verify 
another more recently-proposed dynamic equivalence, this time between the long-time dynamics 
of a Brownian fluid and its corresponding atomic liquid (i.e., the atomic system with the same 
interaction potential) . 

PACS numbers: 23.23.-l-x, 56.65.Dy 
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INTRODUCTION 

At first sight, the macroscopic dynamics of supercooled 
liquids seems to be strongly material-specific, with no 
universal character at all. This is evidenced, for example, 
by the great diversity of molecular glass formers (ionic, 
metallic, organic, polymeric, etc.), giving rise to an over- 
whelmingly rich phenomenology [l|-l5|- One can easily 
understand this lack of universal behavior in terms of the 
wide differences in the materials' structure and composi- 
tion, the masses of their individual atoms, their prepara- 
tion protocol, etc. Of course, the scenario becomes even 
more complex when one attempts to include colloidal sys- 
tems in the discussion. 

The formation of colloidal glasses and gels has been 
the subject of intense study during the last two decades 
0, and it is a widespread notion that the phenomenol- 
ogy of both, the glass transition in "thermally-driven" 
molecular glass formers, and the dynamic arrest tran- 
sition in "density-driven" hard-sphere colloidal systems, 
might share a common underlying universal origin f^. 
Two relevant conceptual issues, however, must be under- 
stood in order for this expectation to have a more fun- 
damental basis. The first one requires us to spell out the 
manner in which undercooling an atomic liquid might be 
equivalent to overcompressing a colloidal liquid. The sec- 
ond is to clarify under what conditions the macroscopic 
dynamics of both classes of systems could be expected 
to be equivalent, given the fact that the microscopic dy- 
namics is Newtonian in atomic liquids and Brownian in 
colloidal fluids. 

The answer to these two questions is highly relevant 
since it will allow us to understand which aspects of the 
macroscopic dynamics of a given system are universal 
and which ones are system-specific. These two issues 
have been addressed using computer simulation methods 



on well defined model systems. For example, interesting 
scalings of the equilibrium dynamics of simple models of 
soft-sphere glass formers have been exposed by system- 
atic computer simulations 0] , which provide an initial 
clue to the possible physical origin of the equivalence be- 
tween the process of cooling and the process of compres- 
sion. Similarly, also using computer simulations, it has 
been partially corroborated that standard molecular dy- 
namics will lead to essentially the same dynamic arrest 
scenario as Brownian dynamics for a given model system 
(i.e., same pair potential) (lol - [l^ . 

From the theoretical side it would be desirable to 
have a unified description of the macroscopic dynam- 
ics of both, colloidal and atomic liquids, which explicitly 
predicts the aspects of the macroscopic dynamics that 
are expected to be universal. These topics might be ad- 
dressed in the framework of a theory such as the mode 
coupling theory of the ideal glass transition [l3|. In fact, 
the similarity of the long-time dynamics of Newtonian 
and Brownian systems in the neighborhood of the glass 
transition, for exam ple, has been studied within this the- 
oretical framework |14|. A number of issues, however, 
still remain open |ll| . 

The present paper is part of an effort aimed at ad- 
dressing these two fundamental issues within a general 
theoretical framework, nam ely, the generalized Langevin 
equation (GLE) formalism [15l4l7|. This formalism was 



employed in the construction of the self-consistent gener- 
alized L ang evin equation (SCGLE) theory of colloid dy- 
namics [18l420l |. eventually applied to the description of 
dynamic arrest phenomena |2ll423l | , and more recently, to 
the construction of a first-principles theory of equilibra- 
tion and aging of colloidal glass- forming liquids [i^, [1^ • 
When applied to model systems with soft repulsive in- 
teractions j26j, the SCGLE theory of colloid dynamics, 
together with the condition of static structural equiv- 
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alence between soft- and hard-sphere systems, predicts 
the existence of a "hard-sphere dynamic universahty 
class" , constituted by the soft-sphere systems whose dy- 
namic parameters, such as the a-relaxation time and 
self-diffusion coefficient, depend on density, tempera- 
ture and softness in a universal scaling fashion 
through an effective hard-sphere diameter determined by 
the Andersen- Weeks-Chandler [H, criterion. These 
predictions provide a more fundamental explanation of 
the scalings previously exhibited by computer simula- 
tions d, Q , and point to the physical basis of the dynamic 
equivalence between cooling and compressing. 

The main purpose of this paper is to report the re- 
sults of, and to provide detailed technical information 
on, a number of simulation experiments performed with 
the purpose of testing this density-temperature-softness 
scaling in the referred dynamic universality class. An il- 
lustrative selection of these results were advanced in a re- 



cent brief communication 27| . The second main purpose 
of the present paper is to perform the pertinent simula- 
tion experiments to test a second relevant prediction of 
the SCGLE theory, which addresses the second of the two 
fundamental issues mentioned above, namely, the macro- 
scopic dynamic equivalence between atomic and colloidal 
liquids. As it happens, the SCGLE theory of colloid dy- 
namics is being extended to describe the dynamics of 
simple atomic liquids 3^, 3l[ . The scenario that emerges 
from these theoretical developments include well defined 
scaling rules that exhibit the equivalence between the dy- 
namics of colloidal fluids and the long-time dynamics of 
atomic liquids. Here we test these scalings by compar- 
ing the simulation results for a given model system using 
both, molecular dynamics and Brownian dynamics sim- 
ulations. 

Thus, the present paper is essentially a report of a set 
of systematic computer simulations, and in Sec. □ we de- 
fine the interaction potentials of the model systems con- 
sidered in our study and provide the basic information 
on the simulation methods employed. In Sec. □ we review 
the concept of static structural equivalence between soft- 
and hard-sphere fluids, and explain how this concept is 
employed to map the static structure of any soft-sphere 
liquid onto the properties of an effective hard-sphere liq- 
uid. In Sec. n we review the extension of this structural 
equivalence to the dynamic domain and present the sim- 
ulation results that validate the accuracy of the result- 
ing dynamic equivalence between soft- and hard-sphere 
liquids. Here we first verify that this dynamic equiva- 
lence is exhibited by the results of our Brownian dynam- 
ics simulations, and then confirm that the same dynamic 
equivalence is also observed in the results of our molec- 
ular dynamics simulations. In Sec. □ we explain the 
correspondence between the dynamics of colloidal fluids 
and the long-time dynamics of atomic liquids, and ver- 
ify that the predicted scalings are indeed satisfied by our 
molecular and Brownian dynamics simulations. At this 



point we have to mention that the present study only 
involves Brownian dynamics simulations that completely 
ignore the effects of hydrodynamic interactions, which 
have an enormous practical relevance in concentrated col- 
loidal fluids. In the last section, besides summarizing the 
main results of this paper, we explain that for the sys- 
tems with interaction potential in the hard sphere dy- 
namic universality class, these effects can be taken into 
account through the value of the short-time self-diffusion 
coefficient, thus expanding the applicability of the scal- 
ings discussed here. 

METHODOLOGICAL ASPECTS 

In this section we describe the most relevant method- 
ological aspects of this work. This includes information 
on the numerical simulation methods and on the theo- 
retical concepts and approaches employed. 

Model potentials 

Let us consider a model liquid formed by N spherical 
particles in a volume V which interact through a soft 
repulsive pair potential u{r) with tunable softness. We 
intend to study the interplay of the effects of the number 
density (or concentration, in the case of colloidal liquids) 
n = N/V, temperature T, and softness, represented by 
some parameter denoted generically as v. There is a va- 
riety of analytic proposals for such tunable soft potential 
, but for concreteness here we shall refer explicitly 
to three specific representative model systems. The first 
is the truncated Lennard- Jones-like (TLJ) potential. 
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in which Q{x) is the unit step function. The positive pa- 
rameter v controls the softness of the interaction, with 
the limit z/ — 00 corresponding to the hard sphere poten- 
tial between particles of diameter a. For fixed v, the state 
space of this system is spanned by the dimensionless tem- 
perature T* = ksT/e and volume fraction <p — 7rncr^/6. 

The specific form of Eq. ([TJ differs qualitatively from 
the inverse power-law (IPL) potential 



u{r) — e{a/ry 



(2) 



commonly used to model hard sphere effects. Just like in 
the case of the TLJ model, for fixed softness parameter ly 
the state space of the IPL model system is also spanned 
by the dimensionless temperature T* = UbT /e and vol- 
ume fraction (f> = Trna^/6. The fundamental difference 
between the IPL potential and the TLJ interaction is 
that the latter is always short-ranged. 
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The third interaction model that we shall refer to is 
defined by the hard-sphere plus repulsive Yukawa (HSY) 
potential, frequently used to model the screened electro- 
static repulsions between charged colloidal particles [s^ . 
This is defined here as 



u(r) 



exp[— z(r/(T — 1)] 



(r/a) 



(3) 



For fixed screening parameter z, the state space of this 
system is also spanned by the volume fraction (j) = 
Trna^/6 and the dimensionless temperature T* = ksT/e 
(sometimes, however, we shall also refer to the repul- 
sion intensity parameter K = l/T* = /3e). The inverse 
screening length z controls the range of the potential, and 
for our purpose, we may consider that it plays the role 
of the softness parameter. Typical values for these pa- 
rameters representing real suspensions of highly charged 
colloidal suspensions at low ionic strength in the dilute 
regime are K = 554, z = 0.149, and (j) of the order of 
10"'* fsS]. We shall use these as illustrative values, along 
with K = 100 and z = 1.0. Figure ([T]) plots these interac- 
tion models for some specific values of these parameters 
to illustrate the variety of interactions considered. 




FIG. 1: Illustration of the TLJ, IPL, and HSY potentials. 



Simulation aspects 

Molecular dynamics (MD) simulations using the 
velocity-verlet algorithm [s^] were conducted for the 
model liquids above, formed by N spherical particles of 
mass M in a simulation box of volume V. The results are 
expressed in the well known Lennard- Jones units, where 
M, a and e are taken as the units of mass, length, and 
energy, respectively, and ^md = Ma"^ / eT* is the cor- 
responding time unit. 

For Brownian dynamics (BD) simulations we follow 
the prescription proposed by Ermak and Mcammon [sS] 
to evolve the positions of the particles in the simulation 



box. Thus, for a given particle at the position r{t) and 
under the force F{t), the particle is displaced in the a 
direction according to 



r„(t + At) = r„(t) + l3D°Fa,{t)At + i?^ 



(4) 



where Z?" is the short-time self-diffusion coefficient of the 
particle. At is the time step, and Ra is a random dis- 
placement extracted from a Gaussian distribution with 
zero mean and variance 2_D"Ai. Taking a as the length 
unit and e as the energy unit, tsD = j becomes the 
natural time scale. 

In both cases, the simulations were conducted with 
N — 1000 particles in a cubic simulation box with peri- 
odic boundary conditions. The initial configurations were 
generated using the following procedure. First, particles 
were placed randomly in the simulation box at the de- 
sired density, such that the maximum overlap between 
particles was in the range 0.65(7 — 0.8ct. To relax this 
initial configuration and reduce or eliminate the overlap 
between the particles we tried two methods. In one of 
them we perform Monte Carlo cycles [35} at a high tem- 
perature, and then decrease the temperature for several 
steps until the original temperature was restored. In the 
other method we uniformly expand the system by in- 
creasing the length of the simulation box by a factor of 
at least 1.5. Then, we run MD or MC cycles while de- 
creasing the simulation box until the original value was 
reached. We checked that the these two methods pro- 
duce equivalent results. Once the initial configuration 
is constructed, several thousand cycles are performed to 
lead the systems to equilibrium, followed by at least two 
million cycles where the data is collected. In the case of 
MD simulations, temperature was kept constant by sim- 
ple rescalingof the velocities of the particles every 100 
time steps [34]. 

Several structural and dynamic quantities are calcu- 
lated from the equilibrium configurations generated in 
the simulations. In particular, the radial distribution 
function g(r) was calculated using the standard approach 
[s^l- The static structure factor S'(fc) can then be ob- 
tained as 



S(k) = 1 -I- 47rn / \g(r) - 1 



sin(fcr) 
kr 



'dr. 



(5) 



Alternatively, S{k) can be calculated directly from the 
positions of the particles in the simulation box 29] . 

Time correlation functions, like the mean squared dis- 
placement (MSD) 



W{t) = {{Ar{t)?)/^ (6) 
and the self-intermediate scattering function Fs{k,t), 
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Fsik,t) 




-ik ■ Arj{t)) 



(7) 



wher 



Af, = 



fj{t) — fj{0), were calculated using the 



efficient, low-memory algorithm proposed in Ref. [36|. 

Crystallinity of the systems was monitored through the 
order parameters Qi, especially Qe, defined as 
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(8) 



where Qim is basically the average, over all particles, 
of the mean spherical harmonics Yim{rij) established 
between each particle i and its close neighbors (j — 
1, iVf,(z)), where Nt{i) is the number of neighbors of 
the particle [s^l- Since in this paper we are interested 
only in the amorphous liquid state, when the simula- 
tions of monodisperse systems exhibited crystalline or- 
der, thus indicating that the corresponding volume frac- 
tion was beyond the freezing point, we discarded that 
monodisperse run, and performed an alternative simula- 
tion introducing size polydispersity to frustrate crystal- 
lization. Polydispersity is handled following a previous 
work [3^, where the diameters of the N particles are 
taken to be evenly distributed between a{l — w/2) and 
(t(1 -I- w/2), with (7 being the mean diameter. We con- 
sider the case w — 0.3, corresponding to a polydispersity 
P = w/^/l2 = 0.0866. Let us emphasize that this pro- 
cedure was followed in both, molecular and Brownian 
dynamics simulations, and that in both cases only size 
polydispersity was introduced, leaving all the other pa- 
rameters unchanged (such as the mass or the short-time 
self-diffusion coefficient of the particles). 

At this point, it is important to emphasize that when 
the system remains in its metastable liquid phase, the 
equilibration time increases enormously as the system ap- 
proaches its dynamic arrest transition (see the detailed 



discussion in Ref. |38|). This means that as the volume 



fraction increases in the metastable region, the initial 
equilibration period will eventually be insufficient, and 
will need to be adjusted to make sure that the system in- 
deed equilibrated properly, as recommended in (38| . The 



present study, however, is not aimed at studying the equi- 
libration process by itself, and hence, we shall avoid ap- 
proaching too close to the glass transition, so as to focus 
our attention on the subject of this work, namely, on the 
dynamic equivalence between soft-sphere liquids. 

To illustrate the end result of this procedure, in Figure 
([2]) we show the corresponding MSD for the TLJ system 
(i/ — 15) in the previous figure for volume fractions in an 
extended range, namely, (j) = 0.1,0.4,0.5,0.6,0.7. These 
results fully cover the stable fiuid phase and much of the 
metastable liquid regime. In all samples the MSD clearly 
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FIG. 2: Mean squared displacement from MD simulations for 
monodisperse TLJ systems with p — 15 and several values of 
the volume fraction (<j!> = 0.1, 0.4, 0.5, 0.6). The most concen- 
trated system is polydisperse, with volume fraction </!> = 0.65. 



exhibits the ballistic and diffusive time regimes typical of 
atomic liquids. We note that the monodisperse system 
with (j) — 0.65 is highly ordered. Hence, we simulate 
instead a polydisperse system whose volume fraction = 
J2iLi ""i^f/S is set to the value cj) — 0.65. 



SOFT-HARD STATIC EQUIVALENCE 

Although simulations are the main methodology em- 
ployed here to generate the static and the dynamic in- 
formation of the model systems above, the analysis of 
this information will rely on a few theoretical notions, 
most notably the predicted static and dynamic equiv- 
alence between soft-sphere and hard-sphere liquids. In 
this analysis, however, we shall recurrently need the ex- 
act structural properties of the fluid of hard-spheres of 
diameter a and volume fraction cj), embodied in its RDF 
9hs{'''/<^]4>) or in its static structure factor SHs{ka; (j)). 
For these structural properties a virtually exact repre- 
sentation is provided by the Percus-Yevick [sol . |40l] ap- 
proximation with its Verlet-Weis correction, defined as 

m 



and 



SHs{kcr;(t>) = S^^^^ka^;^)^), 
with the parameters 0^, and defined as 



0,^ = - </)Vl6, 



(9) 



(10) 



(11) 



and 
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a((/.^/</.)i/3. 



(12) 



The functions g'^^^^ (x; 4>) and S''^^^-' (y; 0) are the so- 
lution of the Ornstcin-Zernike equation with Percus- 
Yevick closure for the HS fluid provided, for example, 
by Wertheim [4^ as easily programmable analytic ex- 
pressions. The resulting 5Hs(r/(T; </)) will be employed 
recurrently in the practical implementation of the con- 
cept of static structural equivalence between soft- and 
hard-sphere systems. This notion was first introduced as 
an essential aspect of the equilibrium perturbation theory 
of liquids [HiiiliH. 

Let us consider a model liquid formed by spherical 
particles interacting through any soft repulsive poten- 
tial u{r) with tunable softness, such as the TLJ po- 
tential in Eq. ([T]). The equilibrium static structure of 
this generic system is represented by the radial distri- 
bution function (RDF) g{r;n,T;a,e,v), also written in 
terms of dimensionless variables as g(r/a; 4>, T* ,1/), with 
T* = ksT/e and cj) = Trna^/6. The physical notion be- 
hind the principle of static equivalence is that at any state 
point (^, T*,z/), this soft-sphere system is structurally 
identical to a hard-sphere system with a state-dependent 
effective hard-sphere diameter ans and effective num- 
ber density nns- This means that for any state point 
{(l),T*,v) one can find a diameter uhs = <^hs{4>jT* ,v) 
and a number density nns = nnsifl^, T*, v) such that 
g{r;n,T;a,e,v) « gHs{r;nHs,o-Hs) H IS Ell , where 
gHs{f",n'HS',o'Hs) is the radial distribution function of 
the HS system, also written as gHs{T/<^HS]4'Hs), with 
4>HS = T^njis^^^sl^- This condition for structural equiv- 
alence can thus be written in terms of dimensionless vari- 
ables as 



{^■,4),T*,i?j « gHs{ 



with 



-^nHSC^HS 



<^HS 



(13) 



(14) 



where A^r is just the state-dependent effective hard-sphere 
diameter in units of ct. 



Xcri(f>,T*,iy) = aHs{4',T*,v)/a. 



(15) 



and A„ is the state-dependent HS particle number density 
in units of n. 



A„(0,T*,i/) = nHs{'i>,T* ,v)/n. 



(16) 



Thus, the condition for structural equivalence can thus 
be written in scaled form as 



g(^^-A,T\u)^gHs(x-'^-:KK<^) ■ (17) 

This equivalence condition can be used in several man- 
ners. The first one is to determine the parameters ans 
and riHS that correspond to a given soft-sphere sys- 
tem at a given state, i.e., to determine the functions 
(THS = cFHs{4>,T* ,v) and urs = nHs{(t>,T* ,v). This 
might be done theoretically, using specific assumptions. 
For example, one could assume that nus = a-nd that 
ghs is 0- independent, with ghs — <^hs{T* ,v) deter- 
mined by means of an approximate version of the equiv- 
alence condition in Eq. (|13p . For example, the approx- 
imation employed in the so-called blip-function method 
reads in general [29| 



47rr 



g-/3tt(r) _ g-^«ffs(i-) 



dr = 0, (18) 



which for the TLJ system can be written as 



\l{T\v)^l-2, I dxx^exp 



(19) 

Evaluating \a{T*,v) determines grs as ansiT* ,1^) — 

Naturally, this or any other approximate scheme has 
a limited range of validity. For example, as we shall 
see shortly, the blip function method is reasonably ac- 
curate for finite-range, moderately soft potentials, such 
as the TLJ liquid with softness parameter 1/ > 6, but it 
fails completely for systems with much softer and longer- 
ranged potentials, such as the HSY fiuid with K — 554, 
z = 0.149, and of the order of 10"''. Thus, it is 
important to search for a more robust method to de- 
termine the functions ans — <^hs {4>i T* , v) and nns ~ 
nHs{4',T* 

One possible method, proposed in Ref. [i^l , is to deter- 
mine the volume fraction of the HS system whose radial 
distribution function provides the best overall fit of the 
exact RDF g{r/a;(j),T*) of the soft-sphere liquid previ- 
ously determined, for example, by computer simulations. 
This method is illustrated here in Fig. ([3]) , where we plot 
simulation data for the RDF of three soft-sphere model 
potentials (TLJ, IPL, and HSY). Fig. for example, 

plots the RDF g{r/a;(j),T*,v) of the TLJ liquid with 
1/ = 6, T* = 1, and (p — 0.7. Thus, we first determine the 
effective hard sphere volume fraction (j^HS by plotting the 
exact RDF gHsir/anS: 4>hs) of Eq. ^ for various vol- 
ume fractions until we identify the value of (j^HS such that 
the height of its second maximum matches the height of 
the second maximum of the soft sphere RDF (w 1.37, 
indicated by the thin horizontal line in the figure). The 
solid curve is the resulting hard sphere RDF. 
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FIG. 3: Radial distribution function for (a) TLJ, (b) IPL, 
and (c) Yukawa systems at the equivalent volume fraction 
(jifig — 0.49. The thick line corresponds to the solution of 
OZ-PY-VW theory for hard spheres. 



This procedure assigns a unique value of (j)HS to that 
set of values of the parameters {(j),T*,i'), i.e., it deter- 
mines the function 4>hs = 4'hs{4'tT* yv). As observed 
in figure (|3^), the height of these two second maxima 
of g{r) coincide, but their positions differ. One finds, 
however, that a simple linear rescaling r — >■ Xa-r of the 
radial coordinate of this HS RDF, prescribed by the 
equivalence condition in Eq. (IT7)) . suffices to match 
the position of the two second maxima. This rescal- 
ing determines the parameter Act, and hence, also the 
effective hard-sphere diameter ans ^t the state point 
{(j),T*,i>) as <7hs = (^Xcr{4>,T* ,v). Finally, the function 
nns = 'nHs(4>i T* , v) is determined by 



riHS 



( 



/ \^HS 



(20) 



Following this procedure in the illustrative example in 
Fig. (|3^) we find that (^hs = 0.49 and \a — (Ths/(^ — 
0.88; as a consequence we find that A„ = nns/n — 1.027. 
These numbers differ only slightly from the results of the 
blip function method, which assumes A„ = 1 and deter- 
mines that Act = 0.888 and <j>HS = 0.49, a comparison 
that illustrates the accuracy of the blip function method 
for the TLJ potential with v — Q. This accuracy improves 
for more rigid potentials and deteriorates for softer and 
longer-ranged ones. 

For example, Fig. ^p) reports an identical exercise 
for the IPL potential with i/ = 3, T* = 1, and </) = 1.2, 
whose RDF is represented by the symbols in the figure. 
Here again the solid line is the RDF of the equivalent 
HS system as a function of t/uhs ^^nd the dashed line 
is the same RDF, but now plotted as a function of r/cr, 
to illustrate the overall agreement between the RDF of 



the soft-sphere system and that of the equivalent HS sys- 
tem. This method determines the effective HS parame- 
ters 4>Hs = 0.49, Act = 0.71, and A„ = 1.007. In con- 
trast, the blip function method, which assumes A„ — 1, 
determines in this case the value Ao- = 1.209 and the 
unphysical HS volume fraction (pus = 2.41. 

Finally, Fig. ([St) reports the same exercise but for 
a much softer and longer-ranged interaction, namely, the 
HSY liquid with K = 554, z = 0.149, and = 2.8 x 10"^. 
As before, the solid line is the RDF of the equivalent HS 
system as a function of r juus- In this case, the resulting 
effective HS parameters are = 0.49, ansl^ — 5.55, 
and TiHs/n = 1.015, leading to the overall agreement 
between RDFs shown in the figure. In this case the 
blip function method (A„ = nns/n = 1) determines the 
completely unphysical values Ao- — <7hs/o' = 27.6 and 
<I>HS = 63.4. Thus, the first conclusion of these three 
illustrative examples is that the assumption that A„ « 1, 
employed in the blip function method above, is indeed 
a very good assumption, at least for these three poten- 
tials at the states structurally corresponding to the HS 
liquid at 4>hs = 0.49. It is then the determination of the 
hard sphere diameter through Eq. (I19p what should be 
avoided, and substituted by the procedure illustrated in 
our examples. 

Let us mention that in each of the three cases corre- 
sponding to panels (a)-(c) of Fig. ([3|) we chose to plot 
the two equivalent RDFs as a function of the radial dis- 
tance r measured in the length unit a of the respective 
system. This comparison, however, can also be done us- 
ing instead the effective HS diameter aus ^ the common 
unit length, as it is done in Fig. ^\X)- There we note, in 
addition, that the simulation data of the three systems 
are actually coincident, and that we only have a single 
HS RDF, corresponding to (pHS = 0.49 and represented 
by the only solid line in the figure. This coincidence illus- 
trates another important feature, namely, that different 
soft-sphere systems that share the same effective HS vol- 
ume fraction also share the same static structure. 

This figure also illustrates the fact that the structural 
equivalence condition in Eq. ()13p can also be used in an 
inverse manner, i.e., to identify the state of a given soft- 
sphere system whose structure matches the structure of 
a prescribed HS system. In reality, what we actually did 
for each of the three soft-sphere model systems in the 
examples in Fig. ([3]) was to search for the state whose 
structure matched the structure of the HS liquid with 
the prescribed volume fraction (j>Hs — 0.49. For this 
we varied the soft-sphere density (or volume fraction (jj), 
keeping the temperature fixed at T* = 1, until meeting 
this condition. 

Another important observation is that for the model 
interaction potentials employed in the present discussion, 
the height of the second maximum of the RDF is not the 
only simple structural order parameter. In reality, some 
other properties that derive from the general equivalence 
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condition in Eq. ([13]) might serve as alternative struc- 
tural order parameters. One of them is the main peak 
of the static structure factor S{k), whose height Smax 
allows us to determine (I)hs , smd whose rescaled position 
kmax determines the effective HS diameter aHS, as illus- 
trated in Fig. dH), which exhibits the structure factors of 
the same systems as in the previous figure, plotted as a 
function, in one case of ka (insets) and in the other case 
of kans (main figure). 



RDF. 




FIG. 4: Structure factors for the TLJ and Yukawa systems 
with effective volume fraction (j>HS ~ 0.49. 

Since our study will extend to densities higher than 
the freezing density of the monodisperse fluid, to study 
the metastable states we need to introduce polydisper- 
sity in our simulation. Although this will not affect dra- 
matically the most relevant dynamic properties, it hap- 
pens to have a profound effect on the thermodynamic 
and structural properties, in particular on the height of 
the second peak of g{r) and of the main peak of S{k). 
These structural order parameters are found to decrease 
with polydispersity (for fixed volume fraction), and this 
requires us to adapt the method described above, to iden- 
tify in a simple manner the effective hard-sphere system 
that corresponds to a given poly disperse soft-sphere liq- 
uid. The adapted procedure is the following. Consider 
a given soft-sphere system with (size) polydispersity P 
and mean diameter a, whose RDF g{r/a] (j),T* ,i>) has 
been measured or simulated. We then model its equiva- 
lent hard-sphere system with the same polydispersity P 
as an equi-molar binary mixture of hard spheres of mean 
HS diameter ans and total volume fraction (I)hs whose 
overall RDF gnsir) (= [ffii + 2gi2 + g22]/2, with gafs be- 
ing the partial radial distribution functions) are obtained 
from the analytic solution of the multicomponent Percus- 
Yevick approximation (43l 44 1 , complemented again with 
the VW correction (i.e., Eqs. (|9))- (IT2]) with a and (p read- 
ing ans and denoted as PY-VW). As in the mono- 
component case, the resulting RDF gHs{f) is then com- 
pared with the simulation results to extract the value of 
ipHS that matches the height of the second peak of the 
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FIG. 5: Simulated radial distribution function g{r) for poly- 
disperse systems with polydispersity P = 0.0866. In the main 
panel we present the results for TLJ systems (symbols) with 
V = & and V = lb, for 4>hs = 0.5. The lines correspond to 
the solutions of the PY-VW theory for the corresponding HS 
systems. Inset: RDFs after scaling the radial distance. 



We illustrate this structural equivalence using a poly- 
disperse version of the TLJ model with different softness, 
u = Q and v = 15, but with densities corresponding to 
the same effective HS volume fraction 4)hs = 0.5, and the 
same polydispersity P = 0.0866 (which is large enough 
to inhibit the crystallization of hard spheres up to very 
high volume fractions (ssjl). Figure ([5]) presents the sim- 
ulation results, along with the theoretical data for the 
corresponding binary mixtures. We see here that the 
scenario depicted in the insets is quite similar to the one 
for monodisperse systems. Furthermore, the scaling of 
the radial distance with the effective HS diameter ren- 
ders the structure of the systems to collapse onto each 
other, as it can be appreciated in the main figure. Thus, 
the results in Figs. ^ and dU show that our proto- 
col to extract the equivalent hard sphere system works 
very well in a wide range of volume fractions for both 
monodisperse and polydisperse systems. With this es- 
sential step covered, we now investigate its implications 
on the dynamics of the equivalent systems. 



SOFT-HARD DYNAMIC EQUIVALENCE 

The dynamic extension of the previous soft-hard static 

Ell 



26 



structural equivalence was discussed in Refs. 
in the context of the dynamics of Brownian liquids, in 
which a short-time self-diffusion coefficient describes 
the diffusive microscopic dynamics of the colloidal parti- 
cles "between collisions". The following discussion also 
refers to Brownian systems, but in the second part of this 
section we shall consider atomic liquids, whose short-time 
dynamics is ballistic. 
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Dynamic equivalence in Brownian liquids 

Although the dynamic equivalence we are about to dis- 
cuss can probably be understood from several perspec- 
tives, in our case we view it as a straightforward pre- 
diction of the self-consistent generalized Langevin equa- 
tion (SCGLE) theory of colloid dynamics and dynamic 
arrest. This theory can be summarized by a closed sys- 
tem of equations for the collective and self intermediate 
scattering functions F{k,t) and Fs{k,t) [2ll - E3l |. which 
in Laplace space read 



F{k,z) 



S{k) 



fc2D0g-l(fc) 

^ ' H-m(fc)AC*(z) 



(21) 



and 



Fs{k,z) = 



1 



l+m{k)AC*{z) 



(22) 



with D'^ being the short-time self-diffusion coefficient. 
These equations become a closed system of equations 
when complemented with the following approximate ex- 
pression for the time-dependent friction function A(^*{t), 



ACit) 



3 {2Try n 



k[S{k) ^ 1] 
S{k) 



dk "^^r/. . F{k,t)Fs{k,t), 

(23) 

and with the following definition of the "interpolating" 
function m{k) [23j 



m(fc) — 



1 



1+1^ 



(24) 



with fi = 2 and with kc being the empirically chosen 
cutoff wave-vector kc = 8.2/a. 

From the condition for structural equivalence in Eq. 
(US]), g{r/a;(f>,T*,iy) « gnsir / cths'Ahs) and assuming 
that A„ = riHsIn = 1 (an excellent assumption in the 
present case), it is not difficult to see that the dimen- 
sionless properties F(fc, i), Fs{k,t), and A(*{z) of a given 
soft-sphere system, can only depend on the wave-vector 
k and the time t through the dimensionless variables 
k* = kans and t* = D^t/aj^g, and that the previous 
equations become identical to those of the hard-sphere 
system at volume fraction (j)HS- This is the origin of 
the predicted dynamic equivalence, summarized by the 
statement that the dynamic properties of the fluid with 
soft repulsive potential u(r), such as the self intermediate 
scattering function (self-ISF) Fs{k,t;n,T;a,e,i^;D°), 
can be approximated by the corresponding prop- 
erty of the (statically) equivalent hard-sphere Brown- 
ian liquid whose particles diffuse with the same 1?°, 



i.e., Fs{k,t;n,T;a,e,i^;D°) w F^"^\k,t;n,aHs; D^). 
This relationship can be written in terms of di- 
mensionless variables as Fs{ka,D^t/ a^\(j),T* ^v) k, 
Ff^\k<jHs.DH/als\'t>Hs).OT 



Fsika, 0, T*,iy) « F^"''\Kka, K^DH/a; A^^), 

(25) 

with Act defined in Eq. (jl5l) . Some consequences of the 
universality summarized by Eq. (j25p were illustrated in 
Refs. [iij and in the context of the TLJ potential. 
Those references, however, discussed in detail only the 
limit of moderate softness {ly ^ 1), in which the strong 
similarity with the HS potential leads to the additional 
simplification that (T/f5(n, T, z^) becomes n-independent, 
and given by the "blip function" approximation fii', ^2^ . 
These, however, are actually unessential restrictions, and 
to illustrate this we have performed Brownian dynamics 
simulations for the non-truncated IPL and HSY models. 

The universality summarized by Eq. (j25p leads to the 
corresponding scaling rules for other properties. For ex- 
ample, let 



Wit;T*,(b,,y) =< (Ar(t))2 > /6 



(26) 



be the mean squared displacement of any soft-sphere liq- 
uid at a given state (T* , 0, v) that structurally maps onto 
the hard-sphere liquid of diameter ans {T* , 0, and vol- 
ume fraction (j)HsiT* , cj), v). Then the normalized MSD 

W*{t*;T*,cj,,v) =<{Av{t*)f > lQals{T\<t^,v), (27) 

with t* = D^t/ajjg, will be identical to that of the equiv- 
alent hard-sphere fluid. 



(28) 



and for that matter, to that of any other soft-sphere liq- 
uid that is structurally equivalent to the HS system with 
the same volume fraction ipHS- 

To test this prediction in Fig. ^ we present the BD 
results for the mean squared displacement of the three 
soft-sphere systems discussed in Fig. ([3]) (i.e., the TLJ 
with v = 6, the IPL with i/ = 3, and the HSY with 
K — 554 and z = 0.149), all of them corresponding to 
an equivalent volume fraction (j)HS ~ 0.49. The MSD is 
presented in the figure in the natural units of the simu- 
lations, i.e. as [W(t)/a'^] vs. [D'^t/a'^]. We observe that 
the MSD exhibits the two linear regimes typical of Brow- 
nian systems 0: at short times [W{t)/a'^] « [D°t/a^] 
whereas at long times [VF(t)/cr^] « D*[D°t/a'^], where 
D* = Dl/D^ is the long-time self-diffusion coefficient 
Dl scaled with Z?°. Thus, at short times the MSD must 
be the same for all systems and states, a condition clearly 
fulfilled by the data plotted in the figure. 
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In the long-time regime, on the other hand, [Ty(i)/cr^] 
is proportional to the scaled long-time self-diffusion co- 
efficient D*, which does depend on the interparticle in- 
teractions, and hence, on the particular system and on 
its state (T*,(/),i/). Thus, the various W{t) in Fig. © 
should differ in their long-time behavior. They, however, 
exhibit the same long-time limit. The reason for this is 
another scaling that derives from the dynamic equiva- 
lence condition in Eq. ((25|). namely, that the dimension- 
less parameter D* depends on (T*, ip, v) only through the 
effective HS volume fraction (pus = (f'HsiT* , (j), v), 



D*{T* A,v) ^ D*hs[(Phs{T* 



(29) 



Since the three systems in this figure were chosen to have 
the same 4)hs {— 0.49), they must share the same value 
of D* (« O.I). According to the scaling indicated in Eq. 
([28| . however, the three systems in Fig. ([3]) must share 
the same W*{t*), i.e., the three different MSDs in the 
figure should collapse onto the same curve when plotted 
as a function of t* = D'^t/aj^g. This is indeed what we 
find, as illustrated in the inset of Fig. which shows 
that the results corresponding to the three curves of the 
main panel do collapse onto each other when plotted in 
this scaled manner. According to Eq. ([28l) . the result- 
ing master curve then determines exactly the function 
W^s[t*'^(t> = 0.49]. The SCGLE theory, besides predict- 
ing this scaling also provides an approximate prediction 
for this function. This prediction is represented by the 
dashed line in the inset of Fig. ([6]), which we see follows 
closely the simulations results. 




We now extend our study to other volume fractions, pay- 
ing special attention to long-time properties such as the 
long-time diffusion coefficient and the a-relaxation 
time Ta- We start by presenting in Fig. ([7^) the BD 
results for D* = Dl/D° (in the format I/D*), as a func- 
tion of (j), for several soft-sphere systems. The main panel 
of the figure, which includes TLJ systems and the solu- 
tion to the SCGLE theory for Brownian systems (i.e., 
Eqs. (I21])-(I24])) apphed to the HS fluids, clearly shows 
how the simulated D* depends on volume fraction and 
softness. For instance, D* is, as expected, a decreasing 
function of 4> and ly. It is also apparent in the figure that 
D* steadily approaches the SCGLE HS results as the po- 
tential becomes stiffer. In particular, all the results in the 
figure converge in the low-(/) limit to the correct limiting 
value D*{(l) ^ 0) = 1. 
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FIG. 7: (a) Scaled long-time diffusion coefficient D* — 
Dl/D^' from DB simulations for TLJ systems with v = 6, 
V = 10 and v — 15, at different volume fractions 0. Inset: 
Results that include the two HSY systems with K — 554, 
z = 0.149, and K = 100 z ^ 1.0. (b) Results plotted as 
a function of the effective volume fraction (f)HS. The solid 
line represents the solution for HS systems from the colloidal 
SCGLE theory. 



FIG. 6: BD results for the mean squared displacement W{t) 
for TLJ {f = 6), IPL {u = 3), and HSY (K = 554, z = 
0.149) systems with the same effective volume fraction (j>HS ~ 
0.49. Inset: Results scaled in HS units. The dashed line 
corresponds to the predictions from the SCGLE theory for 
HS system. 



The results presented in the figure are concerned with 
the full MSD at the single volume fraction ipHS = 0.49. 



Depending on the specifics of the interaction poten- 
tial, other model systems can depart even further from 
the theoretical results. In the inset of Fig. ([7^), for in- 
stance, where we compare the data in the main panel 
with the results for two HSY monodisperse systems, the 
difference with theory and short ranged systems is clear. 
Here the values of D* in the liquid regime are covered 
by the HSY systems for volume fractions in the range 
10^^ — I0~^ (from weakly to highly structured condi- 
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tions), whereas the volume fractions of the TLJ systems 
faU in the typical range 0.1 ^ (/) ^ 1. Despite these differ- 
ences, however, when the data are plotted as suggested by 
the scahng in Eq. (1^ . see Fig. ((TJd), all the data fall on a 
well-defined master curve, which must coincide with the 
exact hard-sphere function _D|jg[0]. In fact, an approxi- 
mation to this function is provided by the SCGLE theory, 
which is represented by the solid line in the figure. Let us 
notice that, although the results in the figure illustrate 
the dynamic equivalence at the very particular condition 
i — >■ 00, basically the same scaling holds at all times 
for time-dependent properties such as the scaled time- 
dependent self-diffusion coefficient D*{t) = W{t)/6D°t 
(illustrated in Fig. (4) of Ref. [il])) or the scaled MSD 
W*{t) (illustrated here in Fig. ([6])). 




FIG. 8: Self-intermediate scattering function Fs{k,t) from 
BD simulations for TLJ systems with = 6 and v = 15, at 
the volume fraction — 0.6. Results are presented for two 
wave vectors fci and fc2, with ki < k2. The solid horizontal 
line represents the target value Fs = 1/e that identifies the 
a-relaxation time Tq of the systems. 

We now turn our attention to the relaxation of the cor- 
relation function Fs{k, t). This is a function that depends 
on the parameters of the system, like the volume fraction 
(/), the wave vector fc, and the softness of the potential. 
In Fig. dH) we illustrate some of its features by plotting 
results for a TLJ system with (p = 0.6, for two values of 
the softness parameter, v = 6 and v = 15, and two wave 
vectors, fcicr = 3.61 and fc2<T = 7.23. The solid horizontal 
line represents the target value Fs{k,Ta) = 1/e, which 
identifies the a-relaxation time Tq, of the systems. The 
data in the figure shows that, for fixed, Fs{k, t) relaxes 
faster with the larger wave vector. The other feature to 
notice is that, for a given value of fc, Fs{k, t) decays faster 
for the softer system (i.e., for v = 6), although the very 
initial evolution is rather similar for the two systems. 
This behavior is expected because the effective excluded 
volume of TLJ systems decreases with the softness of the 
potential, and therefore the structural relaxation is faster 
for softer systems. 

Since the choice for k is not unique, in this work we use 
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FIG. 9: Self-intermediate scattering function Fs{kmax,t) 
from BD simulations for TLJ systems with ly = 6, v = 10 and 
v = 15, and for a Yukawa system with K — 554 and z = 0.149. 
All systems have an effective volume fraction (j}HS = 0.49. In 
the inset the results presented with time scaled according to 
t \-^t. 

k = kfnax^ the position of the main peak of the structure 
factor. In this way structurally equivalent systems will 
then be explored with the same dominant wave vector 
in the hard units, XakmaxO', as indicated by the results 
presented in Fig. ([4]). We evaluate then Fs{kmax,t) for 
two TLJ and one HSY systems with equivalent volume 
fraction (f>HS — 0.49 and the same potential parameter 
CT, and plot the results in Fig. (|9]) as a function of the 
scaled time D^t/a^. In such units, the figure shows that 
the TLJ and HSY systems differ notoriously. In par- 
ticular, the a-relaxation time between the two types of 
systems differs by more than a decade, despite the under- 
lying iso-structurality between the systems, reflecting the 
influence of the potentials on dynamic properties of the 
systems. Dynamic equivalence, as indicated by Eq. (|25l) . 
should however lead to a collapse of all the curves upon 
the transformation to HS units. That this is the case is 
shown in the inset of the figure, where the collapse of the 
curves is achieved by by plotting the data as a function 
of the time scaled not as D'^t/a^, but as D^t/ajjg, as 
suggested by Eq. (|25l). Note here that the other con- 
ditions on the wave vector and the volume fraction im- 
plied by that equation are already fulfilled because all the 
systems share the same effective volume fraction. Such 
scaling of Fs in turn implies that the a-relaxation time 
Tq, a particularly important dynamic property in high 
density or low temperature liquids, turns out to be the 
same for iso-structural systems when expressed in units 
of the equivalent HS system, regardless of the softness of 
the interaction between the particles. 

By performing BD simulations at different volume frac- 
tions we are able to extract the a-relaxation time as a 
function of 0. In the main panel of Fig. ([TO| . where only 
the TLJ systems are presented, Tq. exhibits the typical 
trends found in liquids; it is small for low and moder- 
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FIG. 10: a-relaxation time Tq, in BD units, for TLJ and HSY 
systems at different volume fractions. The figure includes the 
corresponding results for HS systems from the SCGLE theory 
(solid line). Inset: Results that include two HSY systems, one 
with K = 554 and z = 0.149, and the other with K = 100 
and z = 1.0. 

ate volume fractions, and then increases rapidly at a rate 
that depends strongly on the softness of the potential, 
especilly at high volume fractions, as it can be seen in 
the figure. For example, taking as a reference the line 
Ta/tBD ~ Ij the volume fractions reach the approximate 
values (f) = 0.57,0.61,0.77 for v = 15,10,6, respectively. 
Overall, increasing v leads the results progressively close 
to those of HS systems predicted by the SCGLE theory, 
which in the figure are represented by the solid line. 
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FIG. 11: a-relaxation time Ta in HS units. The solid line 
corresponds to the results from the SCGLE theory applied to 
HS fluids. 

In the inset of the figure we replot the results in the 
main panel along with those for two HSY systems. This 
inset then allows us to compare directly the a-relaxation 
time between systems with rather different potentials 
both in range and strength. As expected, the difference 
between HSY and TLJ systems is rather significant, but 
largely dependent on the screening and the strength of 
the HSY potential. Thus, although Tq is haighly sensitive 



to the volume fraction and the conditions of the interac- 
tion potential, the scaling properties of Fs{k,t) already 
illustrated in Fig. ^ should translate to Tq, and make 
the data in the current figure aling with those of equiva- 
lent HS systems. To corroborate this we apply the scaling 
rules for time and volume fraction provided before, i.e. 
we transform the data according to (/) — > (jyns = -^^0 
and [tq]bd = Ta/tBD [toIhs = K'^[Ta]BD, where 
tsD = cr"^ /D^ is the time unit in DB simulations. The 
final results are plotted in Fig. (jlip . where we can ap- 
preciate that the applied transformations work very well, 
collapsing all the data onto a single curve. Note that the 
transformed data follows closely the solid line, i.e., the 
SCGLE theoretical predictions for HS fluids. 

Dynamic equivalence in atomic liquids 

In recent work, our group has extended the colloidal 
SCGLE theory to the domain of atomic liquids. There, it 
has been found that the overdamped colloidal equations 
for F{k,z), Fs{k,z), and AC(0, i-e., Eqs. (EH), 
and p3| are formally equivalent to those of atomic sys- 
tems, subjected to a reinterpretation of the difusion co- 
efficient which in the atomic case depends on density 
and temperature. This formal identity then implies that 
the Newtonian liquids with the same S{k) should exhibit 
soft-hard dynamic equivalence in the same way as Brow- 
nian colloids. In particular, the universality expressed in 
Eqs. (f25|) . (|27| and (|29l) are then extended to molecular 
fluids. Note that such equations are written in terms of 
the scaled time t/tsD, with tsD = o'^/I?°, i.e., in terms 
of the short time property D". Thus, even though the 
same applies to the atomic case, with the correspond- 
ing interpretacion of Z)°, we can choose instead to scale 
the atomic equations with the natural short-time scale of 
Newtonian dynamics, which is Imd = \/ Ma"^ /ksT, thus 
bypassing the explicit form of for atomic liquids. The 
use of Imd as the time scale leads to a different scaling 
rule when transforming to HS units. In this case, time 
is now tranformed as tns = ^athiD, which contrast with 
the Brownian case where tus — X^tBD- 

In what follows, we employ Imd and the correspond- 
ing expresion for tno to study the dynamic equivalence 
between soft-sphere and hard-sphere atomic fluids. It 
should be noted, however, that the static structure of 
our soft systems depends on the volume fraction and the 
speciflcs of the parameters of the potential, but not on 
the underlying dynamics, it being Brownian or Newto- 
nian. As such, the structure factor or the pair corre- 
lation function of a given system, when simulated with 
BD, is identical to that obtained through MD simula- 
tions. In particular, then, the soft-hard static structure 
equivalence stated by equations ([T3| or (fT7|) will still be 
valid for Newtonian fluids, which allows us to discuss the 
soft-hard dynamic equivalence of atomic liquids in the 
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context of iso-structural systems. 



• MD/HS l|)=0.5 
■ MDffLJ v=6 

♦ MD/TLJ v=10 
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FIG. 12: Mean squared displacement from molecular dynam- 
ics simulations for iso-structural soft systems with equivalent 
'f'HS = 0.49. The filled circles correspond to the MD simu- 
lations of a HS fluid. The results are presented in units of a 
and tMD- Inset: MSDs after transforming to HS units. 



share a common long-time diffusion coefficient D l , which 
is defined as the limit DL{(f>) — limt_j.oo W{t)/t. 




FIG. 13: Scaled diffusion coefficient Dl{(j>) = 

D L {<!>)/ {cr^ /tMo) from MD simulations for several soft sys- 
tems. Inset: Results that include the HSY system with 
K = 554 and z — 0.149. The filled circles correspond to 
HS systems. 



We begin by presenting MD results for the MSD for 
some of our working systems, namely, a TLJ system with 
1/ = 6 and 15, a HSY system with K = 554, z = 0.149, 
and an IPL system with = 3, all of them at the effective 
volume fraction c/jhs = 0.49; at the same time, we include 
the corresponding data for the hard spheres fluid, from 
event-driven MD simulations. The MSD for these cases 
are presented in Fig. (fT2)) . where all quantities are scaled 
according to the MD units a and Imd — (Ma'^ /ksT)^/'^ . 
The systems in the figure clearly display the two linear 
regimes, ballistic and diffusive, characteristic of New- 
tonian systems, with the HSY system departing more 
strongly from the exact results (the HS data). Although 
all the systems follow the same curve at short times, the 
HSY system takes longer to cross the ballistic-to-diffusive 
transition, leading to larger values in the MSD at a given 
time. In fact, we have found that our systems follow an 
specific trend depending on the type of potential. We 
have found that, at long times, the MSD of long-ranged 
systems (HSY and IPL) are always above the HS data, 
whereas short-ranged systems (i.e. for TLJ) are always 
below and systematically approaching the exact results 
as the potential becomes stiffer. 

The MSD of these iso-structural systems can be made 
to coincide when plotted in HS units. As before, the 
MSDs are transformed according to equation psp . with 
time scaled in a different manner to reflect the underlying 
Newtonian dynamics of the systems. In the present case, 
the transformation a — > ans leads to t* = t/tMD 
\~^t*, from which we get \^„W{\~^t*). By scaling in 
this manner we appreciate in the inset of the figure that 
all the systems follow nicely the results for the HS system. 
Then we see that all the soft, iso-structural systems have 
the same ballistic and diffusive regimes, and in particular 



The dependence of Dl(4>) on the volume fraction and 
the model interaction is now presented in figure (jl3p in 
the units of the simulation, i.e., as = Dl/ {(j"^ /tMo)- 
The main panel is dedicated to TLJ systems, whereas the 
inset centers on the marked contrast between the strongly 
repulsive, long-ranged HSY system with K = 554 with 
other soft systems, particularly TLJ. Overall, we find 
that the data for mid and high volume fractions exhibits 
similar trends to those already observed in DB simula- 
tions of TLJ systems. The contrast between HSY and the 
other soft systems shown in the inset is also similar. The 
main qualitative difference between MD and DB can be 
observed at low volume fractions, where BD approaches 
the limiting condition £)*(</)-> 0) —> 1 due to the diffu- 
sive nature of its short-time dynamics, whereas from 
MD increases sharply [l/D*^ goes to zero) as the volume 
fraction decreases. In the latter case the Newtonian par- 
ticles perform ballistic motion most of the time, with a 
collision rate that decreases as the system becomes more 
dilute. In fact, the correct </)-dependence of D*^ for dilute 
Newtonian fluids can be obtained from kinetic theory as 
D*]^ ^ l/(t> 49] , which explains the shape of the data at 
low (j). 

Figure (jl4l) summarizes the results for the diffusion co- 
efficient for a variety of soft-repulsive systems (including 
those in the previous figure) , now in units of the equiva- 
lent HS systems. In this case, the applied transformation 
is Dl Df^ = XcDl- Note that we plot l/D^^ and 
to better appreciate that results at both low and 
high volume fractions. From these two forms of the re- 
sults we see clearly that this long-time dynamic property, 
Dl, obtained through the simulation of a variety of re- 
pulsive systems, follows very well the corresponding data 
for truly hard spheres in all the range of volume fractions 
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tentials very close to each other, but far from those cor- 
responding to the HSY potential. More importantly, the 
figure also shows that the Fs{kmax,t) of iso-structural 
systems can be made to fall on top of each other upon 
transforming to HS units. This is illustrated in the inset, 
where we plot the data in the main panel in the scaled 
form Fs{\ak'^g^^, \~^t*), where k'^^^ and t* are the wave 
vector and time in MD units. This result, along with 
the universality expressed by Eq. (|25p . indicates that 
the master curve followed by the systems corresponds to 
the one followed by a Newtonian hard-sphere system at 
the effective volume fraction shared by all the systems, 
(t>Hs = 0.49. 



FIG. 14: Long-time diffusion coefficient DJ, in HS units, as 
a function of the effective HS volume fraction of the systems. 
The left vertical axis presents l/D'^, and the right vertical 
axis DJ,. 



considered in the figure. 



0.6 



0.4 



0.2 











WL » % 0.6 

+1 * J 

i\ » 0.4 






\ \ : 


-f^ — L™ i^naiiiii,. 




^ \ 


10" 10 10" 

a^^(M/kgT)"'"*t 


■ MD/TLJ v=6 

■ MD/TLJ v=10 
A MD/TLJ v=15 
» MD/HSYK=554 

IPL v=3 
— exp 


\ 


V 








10 10 10" 



FIG. 15: Correlation function Fs{kmax,t) from MD simula- 
tions for TLJ [u = 6, 15), IPL {u = 3), and HSY {K = 554) 
systems with the same equivalent volume fraction (j>HS ~ 0.49. 
Time is in units of to = tMD- Inset: Results after trans- 
forming to hard units via to = XatMo ■ The horizontal line 
represents the value Fs{kmax,t = To) = 1/e. 

In what follows we now focus on the properties of the 
atomic self-ISF. As before, the behavior of Fs{k,t) is 
first investigated to show the its principal trends under 
different soft model systems. Later we focus on the long 
time behaviour through the a-relaxation time. Follow- 
ing our previous analysis for Fs{k,t), this function is 
studied at = kmax for several disimilar interaction po- 
tentials. Thus, in figure we present Fs{kmax,t) for 
iso-structural systems, with time and wave vector in MD 
units. The figure includes results for TLJ (i^ = 6 and 
u = 15), IPL {v = 3), and HSY {K = 554) systems, 
all of them with equivalent volume fraction c/jhs — 0.49. 
As we can see, the scenario is rather similar to that of 
Brownian fluids, with the results for TLJ and ILP po- 
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FIG. 16: (a) a-relaxation time as a function of the volume 
fraction, Ta{4i), for selected TLJ, IPL, and HSY systems. The 
filled circles represent the results for HS systems. Inset: Re- 
sults that include two HSY systems, one with K = 554 and 
z = 0.149, and the other with K = 100 and z = 1.0. (b) 
Tc,{(l>) in HS units. 

The above results for Fs{k,t) were presented for iso- 
structural systems at a single value of the equivalent 
hard-sphere volume fraction. Now we address its po- 
tential and volume fraction dependence, as well as to 
its scaling properties, not through Fs{k,t) itself, but 
through its long-time characteristic property Tq . To this 
end, in Fig. ()16b ) we present the a-relaxation time, as 
a function of the volume fraction, for selected soft sys- 
tems. The main panel is dedicated to the softness of the 
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TLJ potential. We observe here that the general trends 
in terms v and </) resemble closely those of Brownian sys- 
tems previusly disscussed. In the inset, where HSY and 
TLJ systems are constrasted, we find that the HSY fiu- 
ids deviate strongly from the TLJ systems, especially for 
the stronger potential (K = 554 and z = 0.149), a similar 
trend found for Brownian systems. The results plotted 
in Fig. (|16b) represent the a-relaxation time for a col- 
lection of soft systems after transforming to HS units, 
i.e., the volume fraction as </) — >■ (t>HS and the relaxation 
time as Ta/tMD — >■ ^aTa/tMD- The results in this figure 
shows that the transformation to hard units works very 
well and that the soft systems follow nicely the MD data 
of strictly hard spheres (black symbols). 



BROWNIAN-ATOMIC DYNAMIC 
EQUIVALENCE 

As mentioned above, our group has extended the SC- 
GLE theory to atomic liquids. There, the overdamped or 
long-time limit of the equations for F{k,t), Fs{k,t), and 
A({t), equations ([21]), and ([23l), have been found to 
be formally equivalent for both colloid and atomic fiuids, 
the fundamental difference being the definition of the dif- 
fusion coefficient D'^. In colloid liquids 1?° is a constant, 
identical to the short-time self-diffusion coefficient given, 
for example, by the EinsteinStokes expression in the ab- 
sence of hydrodynamic interactions. In the atomic case, 
on the other hand, depends on temperature and den- 
sity, and is given by the kinetic-theoretical result ji^ 



Thus, this formal identity implies that the long-time 
dynamic properties of an atomic liquid will then coincide 
with the corresponding properties of a colloidal system 
with the same S{k), provided that the time is scaled as 
_D°i, with the respective meaning and definition of 
One example of such long-time properties is the long- 
time self-diffusion coefficient D^, which in scaled form 
is D* = DL/Do\im^^W{t)/Dot, is predicted to be in- 
distinguishable for atomic and Brownian systems. This 
colloidal-atomic dynamic correspondence was illustrated 
in Fig. (2) of the reference (which we reproduce here in 
Fig. (flT)) ) for some iso-structural systems at the volume 
fractions 4>hs = 0.1 and 0.4. 
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FIG. 18: Long time diffusion coefficient for colloid and 
atomic liquids. In the main panel Dl(</>) is plotted in the 
corresponding units, i.e., in DB or MD units. Inset: Scaled 
results. 
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FIG. 17: Scaled mean squared displacement, simulated with 
molecular dynamics (MD, solid symbols) and Brownian dy- 
namics (BD, empty symbols) for fluids with hard-sphere (HS, 
circles), truncated Lennard- Jones (TLJ, squares) and repul- 
sive Yukawa (Y, trian- gles) interactions, for states corre- 
sponding to (j)HS ~ 0.1 and 0.4. 



This colloidal-atomic dynamic correspondence extends 
over to HS systems and systems with soft repulsive in- 
teractions. We illustrate this in Fig. (|18|) . where we 
plot the long-time diffusion coefficient DL{(j>) for a wide 
range of interacting conditions. The inset of the figure, 
for instance, we plot again our previous results for 
from MD and BD simulations, scaled according to their 
transformation rules. There, one can appreciate the fun- 
damental differences between the two types of dynamics 
(colloid and atomic) when presented in the equivalent 
HS units. The main panel, on the other hand, present 
all those results after scaling the difussion coefficient ac- 
cording to D* = Dl/Dq, where Dq is the usual short 
time diffusion coefficient for DB, but it is the temper- 
ature and density dependent expresion in Eq. pO|) for 
MD simulations. Thus, one can see that combining all 
the scaling rules in this work collapse all the data onto a 
single master curve. 

Finally, in Fig. (2a) of Ref. [sij we showed, for a few 
systems, that the self-ISF follows the colloid-atomic long 
time equivalence. Fig. (2b) of the same reference, on the 
other hand, illustrated the equivalence at the level of the 
a-relaxation time. We provide further data for Tq in our 
Fig. (IT^ . In this figure the inset replots our previous 
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FIG. 19: Inset: a-relaxation time Ta{<j)) for Brownian and 
Newtonian liquids in HS units. Main panel: Data presented 
as r* = k^DoTa. 



data for Tq,, in HS units, for both MD and BD dynamics. 
The main panel, on the other hand, plots (as in Ref. 
(3I1] ) the dimensionless a-relaxation time r* = k^D^Ta 
in HS units, using the corresponding expression for Do 
for MD and DB simulations. This panel shows that all 
the systems are equivalent in the high volume fraction 
regime when plottend in the appropiate scale. 



SUMMARY 

In this work we have presented a simulation study on 
the structure and dynamics equivalence between fluids 
with model repulsive potentials and hard sphere fluids. 
We proposed a method to identify structurally equivalent 
systems that is better than other traditional approaches, 
especially for long ranged potentials. 
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